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ABSTRACT 

Assuming that the observed gamma-ray burst (GRB) rate as a function of redshift is proportional 
to a corrected star formation rate, we derive the empirical distribution of the viewing angles of long 
BATSE GRBs, P om (6>), and the distribution of these bursts in the plane of 9 against redshift, P cm (0, z), 
by using a tight correlation between collimation-corrected gamma-ray energy (E 7 ) and the peak energy 

of vF v spectrum measured in the rest frame (E p ). Our results show that P cm (6) is well fitted by a 
log-normal distribution centering at log6>/rad = —0.76 with a width of <Ji og e — 0.57. We test different 
universal structured jet models by comparing model predictions, P th (8) and P th (9, z), with our empirical 
results. To make the comparisons reasonable, an "effective" threshold, which corresponds to the sample 
selection criteria of the long GRB sample, is used. We find that (1) P th (0) predicted by a power-law jet 
model is well consistent with P cra (9), but P th (9, z) predicted by this model is significantly different from 
P em (6, z); (2) P th (6, z) predicted by a single Gaussian jet model is more consistent with P om ((9, z) than 
that predicted by the power-law jet model, but P th (6) predicted by this model rapidly drops at 6 > 0.3 
rad, which greatly deviates P cra (9); and (3) both P th (8) and P th (9, z) predicted by a two-Gaussian 
jet model are roughly consistent with our empirical results. A brief discussion shows that cosmological 
effect on the E 1 — E p relation does not significantly affect our results, but sample selection effects on 
this relationship might significantly influence our results. 

Subject headings: gamma rays: bursts — gamma rays: observations — methods: statistical 



1. INTRODUCTION 

The phenomenon of gamma-ray burst (GRB) is still of 
a great mystery, although significant progress has been 
made in the recent decade (see reviews by Fishman & 
Meegan 1995; Piran 1999; van Paradijs, Kouveliotou, & 
Wijers 2000; Cheng & Lu 2001; Meszaros 2002; Zhang 
& Meszaros 2004; Piran 2005). It is widely believed to- 
day that the central engines of GRBs power conical ejecta 
(jet) to produce the observed GRBs and their afterglows. 
Sharp breaks and/or quick decays of afterglow light curves 
are regarded as evidence of jets (e.g., Rhoads 1999; Sari, 
Piran, & Halpern 1999). Such evidence is rapidly growing 
up in the recent years (e.g., Bloom, Frail, & Kulkarni 2003 
and the references therein). 

The structure of GRB jet is currently under heavy de- 
bate. Uniform jet model (e.g., Rhoads 1999; Frail et 
al. 2001) and universal structured jet (USJ) model (e.g., 
Meszaros, Rees & Wijers 1998; Dai & Gou 2001; Rossi, 
Lazzati, & Rees 2002; Zhang & Meszaros 2002) are two 
currently competing models. In the framework of the uni- 
form jet model, jet opening angle is assumed to be different 
from burst to burst, but energy distribution within a jet is 
uniform. In the scenario of the structured jet model, it is 
assumed that the energy (and/or bulk Lorentz factor) dis- 
tribution within a jet is a function of the angle measured 
from jet axis. The USJ models suggest that all GRB jets 
have the same geometric structure and the same energy 
distribution within jets. In this model, the observational 
diversity of GRB population is resulted only from differ- 



ent viewing angles. Since there are significant dispersions 
in GRB data, quasi-universal structured jet models are 
also proposed (Zhang ct al. 2004a; Lloyd-Ronning, Dai, & 
Zhang 2004; and Dai & Zhang 2005). Numerical simula- 
tions of jet propagation within collapsar context by Zhang, 
Woosley, & Heger (2004b) showed that the jet structure 
seems to be quite complicated, especially when the part of 
> 0.1 rad is taken into account. Such a jet structure is 
difficult to model by a single component. A jet structure 
with two or more components is likely to be more real. 
The observations of GRB 030329 seem to favor such a jet 
model (Berger et al. 2003). The afterglow lights of this 
burst show two temporal breaks — one occurs at 0.5 day af- 
ter the burst trigger in optical afterglow light curve (Price 
et al. 2003) and the other happens at 9.5 days in the 
radio afterglow light curve (Berger et al. 2003). Berger 
et al. (2003) argued that the interpretation for the two 
breaks requires a two-component jet with 8 Ct i — 0.1 rad 
and f? c .2 = 0.3 rad. Millimeter observations of this burst 
also support this jet model (Sheth et al. 2003). Liang & 
Dai (2004) explained the plausible bimodal distribution of 
the observed peak energy of vF v spectra (E p ) with the sim- 
ilar jet structure model. Huang et al. (2004) and Peng, 
Konigl, & Granot (2004) investigated optical afterglows 
from such a jet model. 

Afterglow observations present more detailed informa- 
tion, but in the late afterglow phase jet structure signa- 
tures may be washed out. Combination of prompt gamma- 
ray emissions and afterglow emissions (especially the early 
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afterglow emissions) may be a powerful way to test differ- 
ent jet models. Zhang ct al. (2004), Lloyd- Roning ct al. 
(2004), and Dai & Zhang (2005) showed that the observed 
E; so — E p relation (Amati et al. 2002; Sakamoto et al. 
2004; Lamb et al. 2005; Liang, Dai & Wu 2004), 9 distri- 
bution, GRB distributions in [9, z)-plane and in (E p , 5 7 )- 
planc, luminosity function, and log N — \ogP distribution 
are roughly consistent with simulation results based on a 
quasi-universal Gaussian-like jet model, where Ei so is the 
equivalent-isotropic energy of GRBs, 5 7 is the observed 
fluence in gamma-ray band, and z is the redshift of GRBs. 
Firmani et al. (2004) constrained the isotropic luminos- 
ity function and formation rate of long GRBs by fitting 
models jointly to both the observed differential peak flux 
and redshift distributions, and found evidence supporting 
a jet model intermediate between universal power-law jet 
model and quasi-universal Gaussian structured model. For 
the uniform jet model their result is compatible with an 
angle distribution between 2° and 15°. The USJ model 
has a power to predict the distribution of viewing angles, 
P th (9). Perna, Sari, & Frail (2003) utilized this power to 
test the USJ models, and found that P th (9) is roughly con- 
sistent with the observed one derived from a sample of 16 
events with 9 known from Bloom et al. (2003). However, 
Nakar, Granot & Guetta (2004) made a further analysis, 
and found that P th (9, z) does not agree with the observed 
result. Liang, Wu & Dai (2004) tested the USJ model by 
simulations based on the £i S o — E p relation and the as- 
sumption of a standard energy budget in GRB jets (Frail 
et al. 2001). In their simulations, they adopted a high 
threshold of gamma-ray fluence, and found that simulated 
P{&) and P(9, z) are consistent with USJ model predic- 
tions. 

The largest GRB sample available so far is the BATSE 
GRB sample. Statistical test with this sample might be 
more reasonable and reliable. More recently, Ghirlanda, 
Ghiscllini, & Lazzati (2004) discovered a very tight cor- 
relation between the gamma-ray energy in GRB jet (E 7 ) 
and E , where E — E p x (l + z) (hereafter GGL-relation). 
It gives an empirical way to estimate the 9 of a burst once 
the E p of the bursts is available. In this analysis, we use 
this relationship to derive the empirical 9 distribution of 
the long BATSE GRB sample, P cm (9), and the distribu- 
tion of these GRBs in the {9, z)-planc, P om (6», z). We test 
different USJ models by comparing the model predictions, 
P th (9) and P th (9, z), with our empirical results. The dif- 
ferences of this analysis from both Perna et al. (2003) and 
Nakar et al. (2004) are summarized as follows. (1) The 
sample we used is the long BATSE GRB (T90 > 2 seconds) 
sample (1213 bursts) but not the current #-known GRB 
sample. The current 0-known GRB sample is too small 
and it suffers greatly observational biases and sample se- 
lection effects (e.g.,' Broom & Lcob 2002; Bloom 2003). 
(2) The observed GRB rate is taken as a corrected star 
formation rate proposed by Bloom (2003). When we com- 
pute the 9 for each burst, we assign a redshift from the 
observed GRB rate model by a simple simulation. This 
GRB rate model is also used to calculate the P{9) and 
P{6, z) predicted by the USJ models. (3) The 9 of each 
burst is calculated by the GGL-relation. Since the E p val- 
ues of most of the bursts in our sample are not available, 
we derive the E p values from a relationship between E p 



and hardness ratios (HR). (4) To ensure that the sample 
should be regarded as a complete one and to make the com- 
parisons between empirical results and model predictions 
reasonable, we use an "effective" threshold according to 
our sample selection criteria. (5)Various USJ models have 
been tested. 

This paper is arranged as follows. In section 2, we 
describe our empirical approach to derive P cm (9) and 
P om (6»,z) for the long BATSE GRBs. In section 3, we 
present the theoretical models of P th (9) and P th (9, z) pre- 
dicted by the USJ models. Results are presented in section 

4, and discussion and conclusions are presented in section 

5. Throughout this paper we adopt Qm = 0.3, SIa = 0.7, 
and H Q = 71 km Mpc" 1 s _1 . 

2. EMPIRICAL APPROACH 
The energy release in the gamma-ray band of a GRB jet 



is 



£ 7 = E iso (l - cos6»), 



and i?i so is calculated by 



iso ~ l + z ' 



(1) 
(2) 



where S 1 is the observed gamma-ray fluence, Dl is the lu- 
minosity distance at redshift z, and k is a factor to correct 
the observed fluence in an instrument band to a standard 
bandpass in the rest-frame (1 — 10 4 kcV in this analysis; 
Bloom et al. 2001). The GGL-relation is, 



E 7i50 - [E p , 2 (l + Z)] 



1.5 



(3) 



where E p ^ = E p /100 keV and E Jt50 = E 7 /10 50 ergs (see 
also, Dai, Liang, & Xu 2004). Although the physics be- 
hind this relationship is still poorly understood, the very 
small dispersion of this relationship makes it a reliable way 
to estimate the value of 9, once the E p and z of a burst 
are well measured. We use this relationship to calculate 
the 9 values of the long BATSE GRBs. Combining Eqs. 
1-3, we have 



9 = arccos[l 



1 £p;t(i 



4tt D\ 2g S 7 , 



(4) 



The E v values of the bursts in our sample are only avail- 
able for some bright GRBs (Band et al. 1993; Ford et al. 
1995; Preece et al. 1998, 2000). It is well known that the 
distribution of E p and HR of long BATSE GRBs arc nar- 
rowly clustered (e.g., Preece et al. 2000). In a previous 
work we showed that the E p and HR of those GRBs with 
a moderate E p (100 - 1000 keV) and HR (1.6 - 6, calcu- 
lated by the fluence measured in the energy band 110-300 
keV to that in 55-110 keV) are strong correlated (Cui, 
Liang, & Lu 2005). Thus, we estimate the E p values by 
the correlation between E p and HR. We show this corre- 
lation for a bright GRB sample of 149 bursts presented by 
Lloyd-Ronning & Ramirez-Ruiz (2002) in the upper panel 
of Figure 1. We perform a linear least square fit at ler 
confidence level to the two quantities. We have 

logEp = (1.86 ±0.14) + (1.16 ±0.09) log HR (5) 

with a linear correlation coefficient of 0.74 and a chance 
probability of p < 10~ 4 . The reduced x 2 is 3.84. We mark 
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the lcr region in Figure 1 (gray band). The distribution of 
the HR for the long BATSE GRB sample is also shown in 
the lower panel of Figure 1. From Figure 1 we find that 
most of the GRBs in our sample well follow this relation- 
ship. 

The rcdshifts of most GRBs in our sample are unknown. 
Please note that our purpose is to examine whether or not 
P em (0) and P em (9, z) are consistent with those predicted 
by the US J models. A GRB rate model has to be used in 
our calculations. To compare the model predictions with 
the empirical results, the same GRB rate model should be 
used to derive the empirical results and model results. We 
assume that the observed GRB rate as a function of red- 
shift for the long GRBs is the same as that derived from 
the currently redshift-known GRB sample. Bloom (2003) 
considered a correction factor of D^ 2 for z > 1 and showed 
that the models of star formation rate as a function of red- 
shift, SF1, SF2, and SF3, from Porciani & Madau (2001), 
are consistent with the observed GRB rate. We hence use 
the same GRB rate model suggested by Bloom (2003) in 
our calculations. Since the largest redshift of GRBs ob- 
served so far is 4.5, we also limit z < 4.5. For a given 
burst, its redshift is assigned by a simple Monte Carlo 
simulation. The procedure of our simulation is as follows: 
(1) obtain the differential probability of the observed red- 
shift distribution with a bin size of 0.01, dQ(z)/dz; (2) de- 
rive the accumulated probability distribution of redshift, 
Q(z) (0 < Q(z) < 1); (3) generate a random number, m 
(0 < to < 1); and (4) assign a z value to the bursts by 
z = (zi+i + Zi)/2, where Q{zi) < m and Q(z i+ i) > m. 
Please note that the redshift assigned by this way for a 
burst can be used only for a statistical purpose. Such a 
simulation ensures that P{&) and P(9, z) predicted by the 
USJ models and derived from BASTE observations are 
based on the same GRB rate model. 

Based on the analysis above, we calculate the 9 values 
by Eq. 4. The uncertainty of 9 is 



de 



de 



" '(1.5^)» 



(6) 

where y = [E Pt2 (l + z )\ /E iso ,50, x = l-y. The uncer- 
tainties of Ei so and E p are given by 



4ttD 2 t 



(7) 



and 



(8) 

respectively, where a = 1.86 ± 0.14 and b = 1.16 ± 0.09. 
Please note when we calculate the 0\E iao , the uncertainty 
of k is ignored since the spectral parameters of the bursts 
are not available (we assume a = — 1 and (3 = —2.3 to 
compute k values for all the bursts). 



3. THEORETICAL MODELS 



The P(6) and P(9, z) predicted by the USJ models for 
a given detection threshold are 



P th (9) =sin<9 



Jo 



dz 



Rgrb dV 
1 + z dz 



(9) 



and 



R, 



""<«•«>— '*nr£. < 10 > 

respectively, where Rgrb is the GRB rate per unit comov- 
ing time per unit comoving volume, dV/dz is the comoving 
volume element at z, and z ma x is the maximum redshift up 
to which a burst with viewing angle 9 can satisfy the detec- 
tion threshold. In this analysis, the Rgrb is taken as the 
observed GRB rate suggested by Bloom (2003). The z max 
is determined by instrument sensitivity and jet structure 
model. It can be derived from 



dV 



/■tt/2 

Jo ' 



4we(9)sm8d8 = L ^ maxJUth (l~cos9), 

1 ~\~ Zmax 

(ii) 

where S t h is the "effective" threshold of the long GRB 
sample observed by BATSE, and e(9) is energy density 
in per solid angle as a function of 9. Please note that 
Sth docs not correspond to the BATSE threshold. An 
observed GRB sample is significantly affected by observa- 
tional biases and sample selection effects, especially when 
the completeness at low fluxes is considered. The "effec- 
tive" threshold should correspond to a given sample selec- 
tion criterion to ensure that the sample should be regarded 
as a complete one in this threshold. Such a sample should 
be selected from a sensitive all-sky survey. The fiuence dis- 
tribution of our sample is shown in Figure 2. It is found 
95% of the bursts satisfy S > 3.2 x 10~ 7 ergs. cm~ 2 . We 
thus take S t h as 3.2 x 10~ 7 ergs. cm~ 2 . The energy den- 
sity profile e(6) for a single power-law jet, Gaussian jet, 
and two-components jet are written as follows. 
Power-law jet: 



^r 2 



(12) 



( wffct ^c is the core energy density when 6 < C . 
^"■©aussian jet: 

e{6) = e e- e2/2e ° (13) 

where 6*o is a characteristic width of the jet, and eo is the 
maximum value of the energy density. 
Two-component jet: 

e(0) = e o (e-* 2 / 2 ^+Ae- e / 2 ^) (14) 

where 6 c ,i and 9 C ^ are respectively characteristic widths of 
the two components, and A is the ratio of energy densities 
in the two components. 
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Based on empirical approach and theoretical models dis- 
cussed above, we calculate the empirical and theoretical 9 
distributions and the distributions of these GRBs in the 
(9, z)-plane. By comparing these empirical results with 
model's predictions, we test these jet models. 

Zhang et al. (2004) showed that the current GRB/XRF 
prompt emission/afterglow data can be described by a 
quasi-Gaussian type structured jet (or similar structure 
jet) with a typical opening angle of 0.1 rad and with a 



4 



standard jet energy of ~ 10 51 ergs. In addition, Liang 
(2004) showed that the fiuence of those GRBs with 9 < 0.1 
rad is almost a constant. We thus take 9q = 0.1 rad and 
e = 10 51 ergs for the Gaussian jet. For the power-law jet 
model, we take 9 C = 0.1 rad and let e c be variable. We find 
that e c = 2.3 x 10 50 ergs yields the best consistency be- 
tween P cm (6>) and P th (9). Bergeretal. (2003) argued that 
the observations of GRB 030329 require a two-component 
explosion with Ci i =0.1 and 9 C , 2 = 0.3. Liang & Dai 
(2004) obtained the similar results based on the E p distri- 
bution observed by BATSE and HETE-2. We thus take 
Ci i =0.1 and 8 Ct 2 = 0.3 for the two-component model. 
The ratio of energy densities in the two components, A, 
is suggested to be 10~ 17 in Liang & Dai (2004). In this 
analysis we let it be variable in the range of (0.01, 0.1) to 
derive a good consistency between P em (0) and P th (#). We 
use different values of A, and find that A <~ 0.08 gives the 
best consistency. The eo in the two-component jet model 
is taken as (10 ± 4) x 10 50 ergs, where the error is assumed 
to be normal distributed. These model parameters are 
summarized in Table 1. 

The empirical results and their comparisons to model's 
predictions are shown in Figs. (3)- (5). The comparisons 
between P em (9) and P th (9) are shown in the left panels 
of the three figures, and comparisons between P em (6*, z) 
and P cm (6*, z) arc shown in the right panels. The error 
bars of P em (9) are calculated by the following method: we 
first derive P em (6» + a g ) and P om (6» -a e ), and then obtain 



CTpcm(fl) by CTpc 



'(<?) 



\P em {9)(9 + a g )-P c ™(9)(9-a 9 )\/2. 



Those GRBs with 0-known from Bloom ct al. (2003) are 
also marked by star symbols in the (6, z)-planes. The dis- 
tribution of these GRBs arc roughly consistent with the 
2(7 region of the P cm (0, z). We quantify the difference be- 
tween P em (9) and P th (9) by a Kolmogorov-Smirnov (K-S) 
test. The result of the K-S test is depicted by a statistic 
of Pk-s- a small value of Pk-s indicates a significant dif- 
ference between the two distributions (Pk-s = 1 means 
two distributions arc identical, and Pk-s < 1-0 x 10~ 4 
suggests that they are significantly different). The K-S 
test results are also listed in Table 1. For the comparison 
between P em (9, z) and P em (9, z), we do not have a quan- 
tifiable way to estimate the difference but evaluate it by 
eye instead. 

From Fig. (3) and Table 1, we find the empirical P om (9) 
is quite consistent with P th (9) predicted by the power-law 
jet model. The K-S test shows Pk-s = 0.575, strongly 
indicating an agreement between P em (9) and P th (9). The 
P em (9) distribution is well fitted by a log-normal function, 



P em {9) = 



0.67 



-2( 



0.57 I 



0.57^ 



(15) 



However, the P cm (9, z) is greatly different from P th (9, z) 
predicted by this jet model. Even in la regions, they are 
quite different. The lack of high-z and large-0 (the right- 
top of the right panel of Fig. 3) make most of this dif- 
ference. We adjust the parameters of this jet model, but 
we still do not get an agreement between P om (#, z) and 
P th (6», z). 

Shown in Fig. (4) is the comparisons of the empirical 
results with whose predicted by the Gaussian jet model. It 
is found that the P th {9) rapidly drops at 6 > 0.3 rad. This 
is significantly different from the empirical result. The K- 



S test for the two distributions shows Pk-s = 2.1 x 10~ 9 , 
indicating that the hypothesis of consistency is rejected. 
Similar to the P th (0, z) predicted by this jet model, the 
P th (9, z) predicted by the Gaussian jet model is also devi- 
ates the P em (#, z). However, comparing the results shown 
in the right panels of Figures 3 and 4, one can find that the 
P th (#, z) predicted by the Gaussian jet model is more con- 
sistent with P em (9, z) than that predicted by the power- 
law jet model. 

The comparisons of P cm (8) and P cm (9, z) with those 
derived from the two-component jet model are shown in 
Figure 5. The K-S test for the distributions of P em (9) and 
P th (9) obtains 3.55 x 10~ 4 , indicating the consistence be- 
tween the two distributions is marginally accepted. The 
P th (9, z) predicted by this jet model is also roughly con- 
sistent with P em (6>, z) in la region. 



5. CONCLUSIONS AND DISCUSSION 

Based on the assumption that the observed GRB rate 
is proportional to a corrected star formation rate, we have 
derived the empirical distributions of long BATSE GRBs, 
P cm (6>) and P om (6>, z), by the GGL-relation. Our results 
show that P cm (9) is well fitted by a log- normal distri- 
bution centering at log# = —0.76 with a width of 0.57. 
We test different USJ models by comparing model pre- 
dictions, P th (9) and P th (#, z), with the empirical results. 
To make the comparison reasonable, an "effective" thresh- 
old, which corresponds to the sample selection criteria of 
the long GRB sample, is used. We find that P th (6>) pre- 
dicted by the power-law jet model is well consistent with 
P om (6>), but P th (9, z) predicted by this model is signifi- 
cantly different from P cm (9, z). Inversely, the P th (9, z) 
predicted by a single Gaussian jet model is more consis- 
tent with P cm (9, z) than that predicted by the power-law 
jet model, but the P th (9) predicted by this model rapidly 
drops at 9 > 0.3 rad, which greatly deviates P em (9). Both 
P th (9) and P th (9, z) predicted by a jet model with two 
Gaussian components are roughly consistent with P cm (9) 
and P cm (9, z). 

The structure of jet is crucial to understanding the na- 
ture of GRBs, such as their rate, luminosity function, and 
explosion energy. Numerical simulations of jet propaga- 
tion within collapsar context by Zhang, Woosley, & Heger 
(2004b) showed that the jet structure seems to be quite 
complicated. A universal, single-component jet seems to 
be difficult to describe such a jet structure, and a jet with 
two (or more) components might be more reasonable. The 
observations of GRB 030329 (Price et al. 2003; Berger et 
al. 2003; Sheth et al. 2003) and the bimodal distribution 
of the observed E p of HETE-2 bursts (Liang & Dai 2004) 
are evidence for this jet structure. The results of this work 
also seem to support this jet model. Zhang et al. (2004a), 
Lloyd-Ronning et al. (2004), and Dai & Zhang (2005) 
showed that quasi-universal Gaussian jet can well interpret 
the observations of the current well follow-up GRB sam- 
ple. Please note that this sample is a bright GRB sample. 
The median of the fiuence in this sample is <~ 2.5 x 10 -5 
erg cm~ 2 (Bloom et al. 2003). It is possible that the emis- 
sions of these bursts are dominated by the core component 
of the jet. Thus, a quasi-universal Gaussian jet might well 
explain the observations of these bursts. 

The true GRB rate as a function of rcdshift is difficult 
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to establish from the current GRB sample. It is generally 
suggested that the GRB rate follows the star formation 
rate. Various models of the star formation rate are pre- 
sented in the literature. Whether or not the model of the 
star formation rate affects our results? We use a model 
of the observed GRB rate suggested by Bloom (2003), 
who constructed a correction factor of when z > 1. 
This factor leads to the observed GRB rate deeply decay 
when z > 1. Thus, different models of the star formation 
rate may give almost the same observed GRB rate (Bloom 
2003), and then the model of the star formation rate does 
not significantly affect our results. In fact, we focus on the 
comparison of both results on theoretical bases and on the 
observational bases. The results of the comparison should 
not be greatly affected by the model of star formation rate 
(Liang, Wu & Dai 2004). 

Our empirical results are derived by the GGL-relation. 
This relationship depends on cosmological parameters. In 
this work we adopt Qm = 0.3 and f^A = 0.7. We 
check if the cosmological parameters significantly affect 
the results of the comparisons between our empirical re- 
sults and model predictions. We take Qm = 0.5 and 
Oa = 0.5. In this case, the GGL-relation becomes E 1 50 = 
(0.90 ± 0.12)£i- 42±01 ° for a GRB sample presented by 
Xu et al. (2005). We show the comparison between the 
P cm (0) based on this relationship and P th (8) predicted 
by the power-law jet model in the cosmology model with 
fl M — 0.5 and f^A = 0.5 in Figure 6. The K-S test for the 
two distributions shows Pk-s = 0.123, confidently sug- 



gesting that they are consistent. Comparing the results 
shown in Figure 6 to that shown in the left panel of Fig- 
ure 3, one can find that cosmological parameters do not 
significantly affect our results. 

We should clarify that our empirical approach and the- 
oretical model are independent without toutology. The 
model predictions are statistical distributions, while the 
empirical results are based on the relationships related to 
the spectral properties and energy release of GRBs. The 
empirical approach and theoretical model are intrinsically 
different. 

Our empirical results are roughly consistent with the re- 
sults from currently 0-known GRB sample. This is a self- 
consistent result because the GGL-relation was discovered 
by this GRB sample. For the bursts in this sample their 
peak energies, temporal breaks of optical afterglow light 
curves, and rcdshifts are well measured. Such a sample 
must suffers greatly observational biases and sample selec- 
tion effects, especially when the completeness at low fluxes 
and the bias of redshift measurement are considered. Band 
& Preece (2005) argued that the GGL-relation may be an 
artifact of the selection effects, and these selection effects 
may favor sub-populations of GRBs. If it is really the case, 
the selection effects should affect our empirical results. 
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Fig. 1. — Upper: log E p as a function of log HR for a sample of 149 bright GRBs presented by Llyod-Ronning & Ramirez-Ruiz (2002). 
The solid line is the best fit line, and the grey band marks the la region of the best fit. Lower: The number distribution of log HR for 1213 
long BATSE GRBs. Two vertical-dotted lines mark a region within which bursts follow the E p — HR relationship shown in the upper panel. 
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Fig. 2.— Number distribution of the observed fluence in energy band 25 - 2000 kcV for 1213 long BATSE GRBs. 
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Fig. 3. — Comparisons between the GRB probability distributions, P(0) (left panel) and P(6, z) (right panel), obtained by our empirical 
approach and by the power-law jet model. In the left panel, the step line with error bars is our empirical results, and the solid line is the 
result of the model. The line contours in the right panel are our empirical P(6, z) (1 a and 2 a) and the grey contours are the results of the 
jet model. The stars in the right panel are the distribution of the bursts with 9 being derived by jet break times (from Bloom ct al. 2003). 
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Fig. 4. — Comparisons between the GRB probability distributions, P(0) (left panel) and P(8, z) (right panel), obtained by our empirical 
approach and by the Gaussian jet model. In the left panel, the step line with error bars is our empirical results, and the solid line is the 
result of the model. The line contours in the right panel are our empirical P(0, z) (1 <r and 2 <r) and the grey contours are the results of the 
jet model. The stars in the right panel are the distribution of the bursts with 9 being derived by jet break times (from Bloom ct al. 2003). 
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Fig. 5. — Comparisons between the GRB probability distributions, P(0) (left panel) and P(8, z) (right panel), obtained by our empirical 
approach and by the two-component jet model. In the left panel, the step line with error bars is our empirical results, and the solid line is 
the result of the model. The line contours in the right panel are our empirical P(6, z) (1 a and 2 a) and the grey contours are the results 
of the jet model. The stars in the right panel are the distribution of the bursts with 8 being derived by jet break times (from Bloom et al. 
2003). 
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Fig. 6. — Comparisons between the 9 distributions derived by our empirical approach (step line with error bars) and by the power-law jet 
model (solid line) in a cosmology with parameters Qm = 0.5 and Qa = 0.5. 



